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Abstract 

Motivated by the TeV-scale gravity scenarios, we study gravitational radiation in the head-on 
colhsion of two black holes in higher dimensional spacetimes using a close-limit approximation. 
We prepare time-symmetric initial data sets for two black holes (the so-called Brill-Lindquist 
initial data) and numerically evolve the spacetime in terms of a gauge invariant formulation for 
the perturbation around the higher-dimensional Schwarzschild black holes. The waveform and 
radiated energy of gravitational waves emitted in the head-on collision are clarified. Also, the 
complex frequencies of fundamental quasinormal modes of higher-dimensional Schwarzschild black 
holes, which have not been accurately derived so far, are determined. 
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I. INTRODUCTION 



Clarifying the nature of black holes in higher-dimensional spacetimes has become an 
important issue since the possibility of the black-hole production in accelerators was pointed 
out. If our space is the 3-brane in large jl| or warped [2j extra dimensions, the Planck 
energy could be of O(TeV) that is accessible with the planned accelerators. If the number 
of dimension D of our spacetime is actually larger than 4, a black hole of very small mass 
will be produced artificially in particle experiments and the evidence may be detected. 

The possible phenomenology of the black holes which may be produced in accelerators 
was first discussed in P| (see ^ for reviews). During the high-energy particle collision 
of sufficiently small impact parameter in a higher-dimensional spacetime, two particles will 
form a distorted black hole of small mass. Subsequently, it settles down to a stationary state 
after emission of gravitational waves. The stationary black hole will soon be evaporated by 
the Hawking radiation, indicating that the quantum gravity effects will be important. The 
evaporation and quantum gravity effects 0, 0] have been studied for yielding a plausible 
scenario (cf. [7] for related issues). On the other hand, the analyses for formation of the 
black hole and for the subsequent evolution by gravitational radiation are still in an early 
stage. These phases are described well in the context of general relativity (see also P|), 
but due to its highly nonlinear nature, the detailed process has not been well understood. 

More specifically, two issues should be clarified for these phases. One is the condition 
(i.e., the impact parameter) for formation of a black hole and the other is the fate of the 
formed black hole after emission of gravitational waves, which can be used as the initial 
condition of the Hawking radiation phase. 

Extensive effort has been made in the past five years for the first issue. The popular 
method is to approximate the high energy particle of no charge or spin by the Aichelburg- 
Sexl shock- wave metric Q|. The merit of this approximation is that superimposing two 
Aichelburg-Sexl metrics, a metric of two particles moving with speed of light can be derived 
for a spacetime region in which causal connection between two particles is absent. Although 
this solution can form a naked singularity at the collision and it is not clear what happens 
after the collision, it is still possible to determine the condition for the formation of an 
apparent horizon for a spacelike hypersurface of the known solution. Formation of the 
apparent horizon is a sufficient condition for formation of a black hole for which the event 
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horizon is outside the apparent horizon. Thus, the lower bound of the impact parameter can 
be estimated. Such study was first done by Eardley and Giddings ^| in a four- dimensional 
case and it was extended to /^-dimensional cases by Yoshino and Nambu . Recently these 
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13| by 



121. 



analyzing the apparent horizon in a 



studies were improved by Yoshino and Rychkov 
different spacelike hypersurface from the one in 

The second issue is to clarify the final state of a black hole formed after collision of two 
particles, which will be a Kerr black hole in higher dimensions, perhaps those described 
by Myers and Perry 3|- Here, the stationary Kerr black holes of no electric charge are 
described by the mass and angular momentum. Namely, the goal is to derive a formula of 
the final mass Mgnai and angular momentum Jfinai as functions of the initial impact parameter 
and initial energy of two particles. For this issue, a couple of preliminary analyses have been 
carried out so far (see Ref. bM for a review). 

n 

One is the work by Yoshino and Rychkov [13||, who constrained the allowed region of 
the mass and angular momentum by finding the apparent horizon for spacetimes of two 
Aichelburg-Sexl particles and subsequently employing the area theorem. However, the al- 
lowed region cannot be pinpointed with this approach, implying that the analysis of gravi- 
tational waves emitted during the collision is inevitable. 

The gravitational radiation from two particles with speed of light in head-on collision was 
first computed by D'Eath and Payne [l6| (summarized in 17 1). They analyzed a spacetime 
of two Aichelburg-Sexl particles in the D = 4 case, paying attention only to a region far 
from the particles and using a perturbative theory. By this analysis, the radiation near 
the symmetric axis can be calculated. Assuming the axisymmetric angular pattern of the 
radiation, they estimated the total radiated energy .Erad as 16% of the total energy of the 
system. 

Recently, Cardoso et al. studied gravitational radiation in the linear perturbation 
theory of the higher-dimensional fiat spacetime. They found again that about 16% of the 
total energy will be emitted in the head-on collision for = 4, which is consistent with the 
results by D'Eath and Payne. They also found that the efficiency is highly suppressed for 
larger value of D, e.g., about 0.001% for D = 10. 

The black hole perturbation theory has been also used recently. Cardoso et al. 
computed gravitational waves from a particle of energy fi with a speed close to the light speed 
falling straightforwardly into a Schwarzschild black hole of mass M ^ fi for D = 4. Berti 
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et al. extend their work for D > 4 (see 



21[ for further generahzations) . Extrapolating 



the resuhs for fi —>■ M/2, they found that the radiation efficiency is ~ 13% for D = 4 and 
decreases with increasing the value of D, e.g., 8% for D = 10. 

Although these approximate studies could give an approximate value of the radiation 
efficiency, it is natural to consider that the error in the estimate is still factor of 2 or more. 
To derive the exact numerical value, it is necessary to carry out more strict analysis. One 
promising approach is to employ numerical simulation in full general relativity. In the D = 4 



case, simulation for black hole collision is feasible [2^, producing certain scientific results. 
However, these works have been done for the case that velocity of each black hole is much 
smaller than the speed of light. Formulation and numerical technique for black-hole collision 
with a very large Lorentz factor 7^1 have not been developed yet. 

In this paper, we adopt the so-called close-limit method for computing gravitational 
radiation, which was originally developed by Price and PuUin 2^. In this method, we 
prepare two black holes of a small separation as the initial condition. If the separation 
is small enough to form a common horizon, the spacetime can be well approximated by a 
perturbed black-hole spacetime. As a result, the gravitational radiation during the collision 
can be analyzed in the context of the black-hole perturbation theory. This method has 



been applied for two black holes initially at rest l2j 

n 
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^initially approaching with linear 
The robustness of this method 



momentum '2d], and many other two-black-hole systems 
is established by confirming that the results by this method agree with those in numerical 
relativity. This fact motivates us to adopt the close-limit approximation for high-velocity 
collision of two black holes in higher dimensional spacetimes. 

As a first step toward the series of study for more plausible cases, in this paper, we 
focus on head-on collision with time-symmetric initial data of two equal-mass black holes. 
For simplicity, we choose the Brill and Lindquist initial data |28] that describes a spacelike 
hypersurface in a spacetime composed of three sheets connected by two Einstein-Rosen 
bridges. We will show that the successful numerical results and indicate that extension of 
the analyses with more general initial data is straightforward. 

This paper is organized as follows. In the Sec. II, we introduce the Brill-Lindquist 
initial data and analyze the apparent horizons for D > 4. In Sec. Ill, we derive the close- 
limit form of the initial data and briefly review the master equation for the perturbation 
of the higher-dimensional Schwarzschild black hole We also explain our numerical 
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methods. Numerical results are shown in Sec. IV, paying attention to the radiated energy 
and gravitational waveforms. Section V is devoted to a summary. In appendix A, the gauge- 
invariant perturbation formalism as well as a method for preparing initial master variable 
from the Brill-Lindquist initial data are presented. Appendix B describes a formula for 
computing the radiated energy of gravitational waves from the master variable. 

II. THE BRILL-LINDQUIST INITIAL DATA 
A. The Brill-Lindquist two-black-hole solution 

Let (S, hab, Kah) denote a (D — l)-dimensional spacelike hypersurface S with the metric 
hab and the extrinsic curvature Kab in a D-dimensional spacetime. The equations of the 
Hamiltonian and momentum constraints are 



where '^""''^^i? is the Ricci scalar of S, is the covariant derivative with respect to hab^ 
and n = D — 2. Assuming the time symmetry (i.e., Kab = 0), Eq. (2) is satisfied triv- 
ially. Assuming further the conformal flatness hab = ^'^^^"'~^^Sab, the Hamiltonian constraint 
equation is written to the Laplace equation for the conformal factor 



where Vf is the fiat-space Laplacian. 

We introduce the cylindrical coordinates (p, z) in which the flat-space metric is given by 
(iSf = dz"^ + dp"^ + p'^dQfj^_i with the metric dQ^^^i on the {n — l)-dimensional unit sphere. 
Among an infinite number of solutions for Eq. ^ that denote spacetimes of two black holes, 
we choose the following one composed of two point sources located at 2; = ±zq along the 
2;-axis as 



- h'^^h'^KacKbd + i^' = 0, 



(1) 



^''{Kab - KbK) = 0, 



(2) 



= 0, 



(3) 
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TABLE I: The critical values of Zq for formation of a common apparent horizon for 4 < D < 11. 



D 


4 


5 


6 


7 


8 


9 


10 


11 




0.192 


0.393 


0.510 


0.586 


0.641 


0.683 


0.715 


0.742 



where R± = =F + p^, M is the total gravitational mass of the system, and r^^M) 
is the gravitational radius defined by 



Here, fi„ = 27r("+i)/Vr((^ + l)/2) is the n-dimensional area of a unit sphere. Hereafter, 
we adopt rh{M) as the unit of the length. The solution ^ provides the system of three 
sheets connected by two Einstein- Rosen bridges. r± = and r oo correspond to spatial 
infinities of each sheet as found by Brill and Lindquist Q for D = 4. 



B. Analysis for apparent horizon 

The close-limit approximation holds for the system sufficiently close to a stationary one- 
black-hole spacetime. Thus, this method can be applied only for the case that a common 
apparent horizon surrounding two black holes is present. Because of this reason, it is nec- 
essary to clarify the range of Zo for which a common apparent horizon exists. Also, by 
obtaining the area of the common apparent horizon, we can estimate the lower bound of the 
final mass as well as the upper bound of the energy radiated away by gravitational waves, 
using the area theorem of black holes. 

The common apparent horizon is determined by a numerical method developed by 
Yoshino and Nambu [s^. Since the system is axisymmetric, it is easily determined by 
a simple shooting method. The shape of the common apparent horizon changes from a 
sphere at = to a spheroid for zq > 0, increasing the ellipticity. At a critical value, ^^o"'*"*; 
it disappears. The values of z^^^^^^ are summarized in Table HI 

In the presence of the common apparent horizon, the mass of the apparent horizon is 
defined by 

^^^=16^ or ' 
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crit) 



FIG. 1: The relation between zq/z^^^^^^ and the trapped energy Mj^n/M at the initial state for 
D = 4-11. The plot range of the vertical axis is 0.9 < Mar/M < 1. 

where Aah denotes the n-dimensional area of the apparent horizon. In the Brill-Lindquist 



initial data, Mah coincides with the Hawking quasilocal mass [31| evaluated on the horizon 
and indicates the trapped energy at the initial state. Mah provides us the lower bound of 
the final mass Mgnai = M — i^j-ad where i^rad is the radiated energy of gravitational waves. 
Equivalently M — Mah gives the upper bound of -Erad- Figure [T] shows the relation between 
2o/4'"*^ and Mah/M. In the D = 4 case, about 99% of the total energy is trapped inside 
the apparent horizon at the initial state. On the other hand, the trapped energy Mah/M 
becomes smaller for larger values of D. 

III. THE CLOSE-LIMIT ANALYSIS 

A. The close-limit of the Brill-Lindquist initial data 

In this paper, the spacetime of two black holes is evolved using the close-limit approx- 
imation, in which the evolution of the field variables are carried out by a gauge-invariant 
perturbation technique. For this analysis, it is necessary to derive an initial condition for 
the linear theory in the Schwarzschild background. To do so, the Brill-Lindquist metric is 
rewritten to 

ds^ = ^^I'^^-^^dR^ + R\de'' + sin^ ednl,^)i (7) 



^ = 1 + 



n— 1 



+ 



— y 



'=2,4,. 



CI 



[{n-l)/2], 



cos( 



(8) 
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where R = a/^;^ + p^, = tan -^{p/z), and Cl^^ denotes the Gegenbauer polynomials which 
are defined by the generating function 



(l-2xt + t2)-' = ^Cfl(a;)t'. 



(9) 



Here, we assume zq <^ rh{M) = 1. If zq = 0, the metric provides the space component of 
the Schwarzschild metric in the isotropic coordinate. 

We introduce a new coordinate r through the ordinary relation between the Schwarzschild 
coordinate r and the isotropic coordinate R: 



r = R^, 



2/(n~l) 
' 



% = l + l/AR 



n-1 



or equivalently 



R 



.(n-l)/2 + V^n-1 _ A 12 



2/(n-l) 



Then the metric becomes 



2 , 



dr 



(10) 



(11) 



(12) 



^)'c;("-^)/^](cos0). 



(13) 



Here, the metric in the square brackets of Eq. ()12|) denotes the space part of the 
Schwarzschild metric. The difference of (P'/i^o from unity is of O(e^) where e = ZQ/R. Since 
the region of r > 1 corresponds to R > R^ = 4-^/*^"-^^ the system can be regarded as the 
Schwarzschild black hole plus its perturbation for a sufficiently small value of zq (or e). 

The first order perturbation includes the mode I = 2, 4, ■ ■ ■ whose order is O(e'). We only 
consider the leading O(e^) correction which is the / = 2 mode. Then we find the prefactor 
(i^/!Z/o)4/(n-i) of Eq. (O becomes 

4/(n-l) 

(14) 



Wo 



1 I l/{n~l)R-' /^o^|V[(n-i)/2], 
l + l/4i?"-i \r) ^ 



cos I 



This provides the major parts of the initial data for the linear perturbation theory in the 
close-limit approximation. 
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B. Master equations and initial master variables 



In order to analyze the time evolution of the perturbation, we use a gauge invariant 
formulation of the higher- dimensional Schwarzschild perturbation 2^. The master equation 
for the three types of perturbation variables, i.e., the scalar, vector, tensor variables was 
derived in Ref . fZ^ . Since the Brill-Lindquist initial data is axisymmetric with no rotation, 
we only need to evolve one master variable of the scalar mode. 

In the scalar-mode perturbation, the master variable $ that is related to the gauge- 
invariant quantities obeys the master equation 

^2$ ^2^ 



where 



dt'^ drl 



+ Vs^ = 0, 



Vs{r) 



f{r)Q{r) 



(15) 



(16) 



/(r) = 1 — X, H{r) = m+ (l/2)n(n + l)x, 



(17) 



m = k'^ — n, k'^ = l{l + n — 1), x = 1/r 



n-l 



(18) 



g(r) = n'^{n + 1) V + n{n + 1) [A{2n^ - 3n + A)m + n{n - 2)(n - A){n + 1)] x^ 

- I2n [{n - 4)m + n{n + l){n- 2)] mx + 16m^ + An{n + 2)m^. (19) 



denotes the tortoise coordinate defined by 



r* = / dr/ f. 



More explicitly. 



(20) 



= r 



n 



n/2-l 

2 X ^ zmn 

> sm 

- 1 ^ n-l 

m=l 



2m7r 2m7r 
arctan | — cot r + r esc 



n — 1 



n — 1 



7r/2 



+ 



1 



n — 1 



n/2-l 

log(r — 1) + cos 



m=l 



2m7r / 2 2m7r 

log 1 + r — 2r cos 

n — 1 V n — 1 



(21) 
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FIG. 2: The initial condition for the master variable ^ jz^ as a function of r,,. The cases = 4, 6, 
8, and 10 are shown. 



for even n and 



r* = r 



(n-3)/2 



2rm\ 



sm ■ 



+ 



m=l 
1 



n 



n 



log 



- 1 

r — 1 



arctan — cot 



2m7T 



n 



r + 



(n-3)/2 

I ) + E 



r?i=l 



1 

2m7r 
— 1 



2m7r , 

r CSC ] — it/2 

n — 1 



loff ( 1 + r 



2r cos ■ 



2m7r 
n — 1 



(22) 



for odd n. 

The initial condition for the master variable $ in gauge invariant perturbation theory 
is calculated from the metric ()12|l and ()14j) numerically (see Appendix A for details). The 
initial values of ^/zq for D = A, 6, 8, and 10 are shown in Fig. |21 ^/zq asymptotes to a 
nonzero value for r^, — — oo and to zero for r^, — oo. It rapidly changes around ~ 0. 



C. Numerical methods and numerical error 



In the numerical computation for Eq. (15), we used the Lax-Wendroff scheme which is 
second-order accurate in time and space. To validate our code, the following convergence 
test was carried out: We evolved the Gaussian wave packet using the initial condition as 
/i(0,r*) = exp(— r^/100) with various grid resolution. Denoting the grid spacings of r* and 
t by dx and dt, respectively, there are four free parameters in our numerical code: dx, 
X = dt/dx, and the locations of the inner and outer boundaries ri™-* and ri"^^\ First we 
computed the fiducial solution hj\f = ^/zq {N denotes the step size of the t direction) at 
r* = 100 for < t < 250 in the D = 4 case, choosing dx = 10'^, X = 0.2, r^"^^ = -200, 
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FIG. 3: The relation between the grid spacing dx of and the numerical error. Our code shows 
the second-order convergence. 

and ri°"*^ = 200. Then, we repeated the computation choosing the larger values of dx while 
fixing the other parameters, and estimated the error by 

Err = ^N\hN-h^\ ^^s) 

for each value of dx. Figure El shows the relation between log^g dx and log^Q Err. All points 
are located on a straight line of which slope is two. This illustrates the second-order accuracy 
of our code. 

Then we analyzed the time evolution of the Brill-Lindquist initial data by solving the 
master variable ^/zq starting with the initial values shown in Fig. |21for A < D < 11. We 
used dx = 0.01 for D = 4-7 and dx = 0.005 for D = 8-11. For the larger values of D, we 
choose the smaller grid spacing since the error increases with increasing D. We chose the 
other parameters to be A = 0.2, rl""*"* = 1000, and ri™^ = —200. Comparing the results with 
those computed in poorer resolutions, we found that the numerical error is within 0.02% for 
all values of D. 



IV. NUMERICAL RESULTS 

A. Waveform and radiated energy 

Figure m shows the behavior of ^/zq for D = 4-11 observed at = 100. The quasi- 
normal mode ringing is seen for all values of D. The power-law tail is also computed well 
for D = 4 and odd D. The reason that the tail cannot be seen for even values of D > 6 is 
that the tail decays more rapidly than that of odd values of D, as clarified in js^. 

11 




The total radiated energy E^ad is calculated by 

E^_^ k\n-l){k^-n) f 
M ^ 2n2n„ J ^ '''' 

(see Appendix B for a sketch of the derivation). Since $ is proportional to for the I 
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TABLE II: The values of E2 = E,^a/4 and E^rad at z = for L> = 4-11. 



D 


4 


5 


6 


7 


8 


9 


10 


11 


i?rad(4'"*^) (%) 


0.0252 
0.0034 


0.0245 
0.059 


0.0290 
0.20 


0.0288 
0.34 


0.0258 
0.44 


0.0224 
0.49 


0.0195 
0.51 


0.0172 
0.52 




FIG. 5: The energy spectrum of gravitational waves. The unit of the vertical axis is Mzq. The 
location of the peak shifts to the right-hand side as the value of D increases. 



mode, i^rad is written as 



E,ad/M^4(^0)'. 



(25) 



Table ITTl shows the values of E2- In the D = A case, E2 has been already obtained by 



Abrahams and Price 



24| as E2 = 0.0251. This agrees well with our numerical result. 



To compare the radiation efficiency Ej-g^^/M among the different values of D, one has 
to specify the values of zq. In Table m we summarize the values at the critical values 



(crit) 



for formation of the common apparent horizon. It is found that i?rad(-2o"'*^) 



increases by increasing the value of D. However we also should mention that the higher- 



order correction might be large for Zq 



(crit) 



. As we can see from Eq. ()14|) . the characteristic 



value of the first order perturbation is {^/^q)^^^'^~^^ — 1, which becomes maximal at the pole 
on the horizon. Such a maximal value is quite large, e.g., ~ 1 for D = 4 and ~ 6 for 
D = 10. Although the close-limit method gives a fairly good approximation beyond the 
regime of the perturbation in the four- dimensional case [2^, further investigations such as 
the second-order analysis or the full numerical simulation are necessary to clarify this point 
in higher- dimensional cases. 
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FIG. 6: The angular dependence A(0) of the radiated energy. There is a peak at 9 = 7r/2 and 
A(7r/2) becomes larger as D increases. 



TABLE III: The values of evaluated at zq < 1 for -D = 4, 11. 



D 


4 


5 


6 


7 


8 


9 


10 


11 


an 


0.0034 


0.016 


0.024 


0.024 


0.021 


0.017 


0.014 


0.011 



Figure El shows the energy spectrum of gravitational waves. The value of uj at the peak 
becomes larger as the value of D increases. This reflects the fact that the real part of 
the fundamental quasinormal mode frequency increases with D (see Sec. IV C for the 
quasinormal mode frequencies). 

The angular dependence of the radiated energy is given by 

E dd ^ ' v^r(2 + n/2) ' ^ ^ 

(see Appendix B for a derivation). Figure IHl shows the behavior of the function A{6) for 
D = 4, 11. Gravitational waves are mainly emitted around the equatorial plane and this 
tendency is enhanced for larger D. This reflects the fact that there are more directions 
transverse to the symmetry axis for larger values of D. 

B. Relation between E'^ad and Mah 

As already mentioned, M — Mah provides the upper bound of -Erad- Hence the ratio of 
E,^ to M - Mah 
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FIG. 7: The relation between M - Mar and ^rad for D = 4,..., 7 (left) and for D = 7, ...,11 
(right). The values of £'rad/(-^ — -^ah) are much smaller than unity. This is consistent with the 
area theorem. 

should be smaller than unity and it provides one consistency check of our calculation. For 
small value of zq, both Ej-ad and M — Mar are proportional to Zq and an take nonzero values. 
The values ofao are summarized in Table ITTTl and then we confirm that an is less than unity. 
The values of ao in higher dimensions are larger compared to for four dimensions. The 
relation between M — Mar and i^rad for < 2:0 < 2:0'^"*'* is also depicted in Fig [3 <C 1 
also holds for this parameter range. 

C. Quasinormal modes 

From the ring down phase seen in Fig. |31 it is possible to derive the complex frequencies 
of the fundamental quasinormal modes cuqn- By comparing them with previous studies, we 
check the reliability of part of our results. In addition, the values of cjqn that have not been 
accurately computed for large values of D so far are derived from our numerical results. 

The values of wqn are evaluated in the following manner: The imaginary part Im(u;QN) 
is derived from the slope of the peaks of log $(t) shown in Fig. 4. The real part is estimated 
by averaging the intervals of zeros of $(t) and consistency is checked by identifying the 
Fourier peak of $(t) x exp(— Im(co'QN)t). For D = 9, 10, and 11, we found that two modes 
are mixed and searched two values of cuqn so that the numerical data of $(t) is well fitted. 
By comparing the results derived by using the several ranges of t, we estimate the error to 
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TABLE IV: The fundamental quasinormal mode of the scalar gravitational perturbation in the 
Schwarzschild black hole for / = 2 and 4 < D < 11. Results by our method, by Leaver's method, 
and by the WKB method are shown. 



D 
4 
5 
6 
7 



Our estimate 



Leaver's method 



WKB 



10 



11 



0.747 - 0.177i 
0.947 - 0.256i 
1.139 -0.305i 
1.339 - OAOOi 
1.537 - 0.587i 
1.19-0.95i 
1.98 - 0.90i 
1.25 -0.94i 
2.47 - 0.99i 
1.20 -0.98i 
2.91 - l.lli 



0.7473 - 0.1779i 
0.9477 - 0.256H 



0.746 - 0.178i 



1.131 - 0.386Z 



[1.778 - 0.571i] 



[2.513 - 0.744i] 



£1=9 




D=10 




-1 

-2 
-3 
-4 

-6 
-7 



£1=11 




104 106 108 110 112 114 
t 



105 110 115 120 125 
t 



102 104 106 108 110 112 114 116 
t 



FIG. 8: The behavior of logj^o ^{'^)/^o ringing phase (the black line) and the fitted data with 

two values of wqn listed in Table HVI fthe gray line) for D = 9, 10, and 11. 



be < 1% for 4 < D < 8 and ~ 5% for 9 < D < 11. We summarize the values of wqn for 
D = 4-11 in Table |TV| and compare the fitted data and $(t) for D = 9, 10, and 11 in Fig. |H1 
The par t of the derived results can be compared with previous ones. For = 4, 
Leaver 33 1 gives very accurate values of cjqn (pe sJj for a review). For the higher- 
dimensional Schwarzschild black holes, Konoplya |35|] evaluated the values of cuqn for / = 2 
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TABLE V: The same as Table IV but for / = 4. 



D 


Our estimate 


Leaver's method 


WKB 


4 


1.618 - 0.188i 


1.6184 - 0.1883Z 


1.618 - 0.188i 


5 


2.193 - 0.328i 


2.1924 - 0.3293i 




6 


2.623 - 0.439z 




2.622 - 0.438i 


7 


3.012 - 0.534i 






8 


3.389 - 0.631i 




3.401 - 0.645i 


9 


3.779 - 0.734i 






10 


4.176 -0.838i 




4.223 - 0.841i 


11 


4.595 - 0.950i 







TABLE VL The same as Table IV but for / = 6. 



D 


Our estimate 


Leaver's method 


WKB 


4 


2.455 - 0.192i 




2.424 -0.191i 


5 


3.286 - 0.344i 






6 


3.913 - 0.470i 




3.911 - 0.467i 


7 


4.610 - 0.574i 






8 


4.924 - 0.674i 




4.923 - 0.675i 


9 


5.388 - 0.769i 






10 


5.834 - 0.859i 




5.848 - 0.865Z 


11 


6.292 - 0.955i 







and 3 using the WKB method and they were extended to / > 4 in 20 1 . The Leaver's method 

" n 

was apphed to the scalar mode of the gravitational perturbation for D = 5 jSy]. The values 
of wqn derived with these methods are summarized in Table HVl In the cases D = 4 and 5, 
our results agree well with those in the previous studies. On the other hand, for D = 6, 8, 
and 10, our results disagree with the previous ones. As stated in HQ, the WKB method 
is not expected to work well for D = 8 and 10 since the potential Vs has a negative peak. 
Our results indicate that the WKB method might not be good even for D = 6. 

Our methods can be applied to arbitrary values of I. To estimate cuqn, we evolved 
appropriate initial data for / = 4 and 6. For these cases, only one quasinormal frequency 
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appears in the ring down phase and the error of cuqn is < 1% for all values of D. Our results 
together with those in the previous studies are summarized in Tables fVl and IVll They agree 
well within the difference < 1%. This implies that evolving appropriate initial data by the 
master equation is an effective method for computation of the complex frequencies of the 
fundamental quasinormal modes. 

V. SUMMARY 

In this paper, we have studied gravitational waves emitted during head-on collision of two 
black holes in higher dimensions using the close-limit approximation. We evolved the Brill- 
Lindquist initial data perturbatively, using a gauge-invariant techniques in a Schwarzschild 
black hole and calculated the waveform and radiated energy -Erad- -Erad is given by the 
formula (|25|) and the values of E2 are summarized in Table IHI At the critical separation 
for the presence of the common apparent horizon zq = zl)^^^^\ our analysis of the first order 
perturbation predicts that Ej-ad/M becomes larger with larger values of D. There is a 
possibility that the higher-order correction is large for zq = ^o"'*"* in the higher-dimensional 
cases and to clarify the higher-order effect is left as a remaining problem. We also evaluated 
the values of an = E^gj^/i^M — Mah) at 2:0 <^ 1 and found that an {5 < D < 11) is 
larger than 04. These results indicate that more energy could be radiated away in higher 
dimensional spacetimes than in the four dimensional one during head-on collision with the 
approaching velocity much smaller than the speed of light. 

It has also been illustrated that the fundamental quasinormal mode frequencies for the 
scalar-mode perturbation in the Schwarzschild black holes can be computed in our analysis. 
We derived the complex frequencies for various values of D, which have not been accurately 
computed so far. 

As we mentioned in the introduction, the close-limit analysis for head-on collision of the 
black holes performed in this paper is the first step toward more general studies of the black 
hole collision in higher dimensions. As the next step, we plan to analyze the evolution of two 
initially-moving black holes using the close-slow approximation. Both the head-on collisions 
and the grazing collisions should be studied. As the momentum of each black hole increases, 
dependence of the emissivity of gravitational waves on the value of D would be changed 
from the results in this paper. By observing such behaviors, we will be able to discuss the 
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dependence of gravitational radiation in the collision process on the value of D in a different 
way from the previous studies. 
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APPENDIX A: DETERMINING THE INITIAL MASTER VARIABLE 

In this section, we explain how to match the initial master variable $ to the initial data. 
We begin by briefly reviewing the gauge invariant formulation of the perturbation in the 



Schwarzschild black hole spacetime 



29|. The background metric is given by 



ds' = gabdy^dy" + r\y)-iijdz'dz\ (Al) 

where Qahdy'^dy^ = —fdt^ + f^^dr"^ and 'jijdz'^dz^ denotes the metric on a unit sphere. The 
perturbed metric is written as 

ds"^ = {gab + hab)dy''dy^ + {hai + h,a)dy''dz' + {r^j + hij)dz'dzK (A2) 

The perturbation variables hab-, hai, and hij can be separated into scalar, vector, and tensor 
modes. Each mode is further expanded into the modes of different anuglar quantum number 
/ = 0, 1, ■ ■ ■ using the hyper-spherical harmonics §, which is the solution of the following 
equation: 

{D^D' + A;2)§ = 0. (A3) 

Here denotes the covariant derivative on the unit sphere and the definition of k"^ is given 
in Eq. ()18|). The variables of the scalar mode perturbation are given by 

hab = fab^, (A4) 

hai = rfaSi, (A5) 
hij = 2r^ (HLiijS + HrSi, ) , (A6) 
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where 



1 



1 



n 



In the axisymmetric case the metric of a unit sphere is written as 



and then 



jijdz'dz^ = de'^ + sin^ 9dQ,l_i, 



§ = 5fi = ii:|"icK"-')/'](cos^), 



47r("+i)/2r(n + 1-1) 

(n + 2l- l)r(l + l)r((n - l)/2)r(n - 1) 



-1/2 



Ki is the normahzation factor that is chosen so that 



J — Sill 



is satisfied. For the gauge transformation generated by the following vector fields, 



the gauge invariant quantities of the perturbation are given by 



F = Hl + {l/n)HT + (l/r)DVX„, 

Fab — fab + DaX^, + D^Xa, 



where 



The master variable $ is related to the gauge invariant quantities as follows: 



X =r"-2(F/-2F) 



/ 



y = r"-2 (F; - 2F) = r"/2-2 (^jdt^ 



16H^ AH 
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where 



Px{r) = n^{n + ifx^ + 2n{n + l)[2(n^ + n + 2)m - n{n - 2){n + l)]x^ 

-An[{n - ll)m + n{n + l){n - 3)]mx + 16m^ + 8m^n^, (A20) 

Qx{r) = n{n + 1) V + 2[(3n - l)m - n{n + - 4nm, (A21) 

Py[r) = n^{n - l){n + ifx^ + 2n{n^ - l)[4m - n{n - 2){n + 

+4ra(n - 1) [3m + n{n + l)]mx, (A22) 
l){n + l)x^ - 2{n - l)[m + n{n + l)]a;, (A23) 

Pzij') = [—n^iji + l)x + 2(n — 2)m]y + n{n + 

+ [2(2n - l)m + n(n + l)(n - 2)]x - 2nm. (A24) 

Here, the definition of x is given in Eq. (fTHjl . 

Now we turn our attention to the method for determining the initial master variable $ 
from a full nonlinear time symmetric initial data. We denote fab and fa as 

/a6=|^^° |, rfa = {hoM- (A25) 

\^ f-^H, ) 

Because of the time symmetry of the initial condition, 

Ho = H2 = h = Hl = Ht = 0, (A26) 

Hi = ho = 0. (A27) 
By the comparison with Eqs. (fT^ and (fT^ . we find 

h^ = Ht = 0. (A29) 
Then the gauge-invariant quantity is found to be 

F = Hl, Fl = 0, (A30) 

and thus 

X + Y = -nr''-'^x{r), (A31) 
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where we have used one of Einstein's equation = — 2(n — 2)F. On the other hand, X + Y 
is given in terms of $ as follows: 



Hence we find the following equation for the initial master variable $ 



/ 



Px + Py 
AH 



dr^ r{Qx + Qy] 

We also find $ = from the condition Z = 0. 

Taking the limit r^, ^ — oo of Eq. ()A33|) . we find 



2 . 4^/(^-^)nzh 



(A32) 



(A33) 



(A34) 

{n-l){n + m)KP 

This gives us the boundary condition at r^: = — oo. We solve Eq. ()A33|) using the fourth- 
order Runge-Kutta method with the boundary condition (A34). 



APPENDIX B: FORMULA OF RADIATED ENERGY 

In this section, we sketch the derivation of Eq. ^M . The similar calculation for the 

n 

energy spectrum is found in Ref. j20||. 

In the region far from the source, outgoing waves are well approximated by the spherical 
ones in the transverse-traceless (TT). Then, the perturbation is written as 

hj;" ^ 2r'HTSij, (Bl) 



HT-^^,hit-r). (B2) 



The radiated energy flux is given by 



^ /.J.-/.™^- = -^ifX-S-, (B3) 



dSdt 327iG' 87cG 
which is the same form as in the four- dimensional case. Using the formula 



we find 



dE Hl {n-l){e-n) 
dt ~ SnG nk^ ■ 
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Next we rewrite the formula ()B5|) in terms of the master variable $. In a region far from 
the source, the gauge invariant quantities become 

F = -Ht + -Xr, Fab = DaXb + D^Xa, (B6) 

n r 

Xa = ^daHT. (B7) 

If we calculate F, F/, F^, and F[ keeping only the leading order 0(r^~"/^) and the subleading 
order 0(r^~"/^), we find 

Y + Z = n-g^Hr, (B8) 

where Eq. flB2|) was also used. On the other hand, calculating Y + Z in terms of $ and 
using the fact that $ = — holds for the outgoing wave, we obtain 

Ht = — r-"/2$. (B9) 
2 



Substituting this equation into Eq. ()B5|) . we find 

^e(n-l)(k'^ -n) /••2, , 

This formula is equivalent to Eq. ()24|) in the unit rh{M) = 1. 
Using Eqs. (jHSI), and (IHTn|l . we find 

In the case / = 2 it becomes 

^j£^2.-(.«,/^r((„ + 5)/2) 

Edn„ n{n + 2) ' ^ ^ 

which reduces to Eq. using = Qn-i sin"~^ 9d9. 
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